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Abstract 

We study 3d CFTs with an 0{N) global symmetry using the conformal bootstrap 
for a system of mixed correlators. Specifically, we consider all nonvanishing scalar 
four-point functions containing the lowest dimension 0{N) vector (pi and the lowest 
dimension 0{N) singlet s, assumed to be the only relevant operators in their symmetry 
representations. The constraints of crossing symmetry and unitarity for these four- 
point functions force the scaling dimensions (A(^,A<j) to lie inside small islands. We 
also make rigorous determinations of current two-point functions in the 0(2) and 0(3) 
models, with applications to transport in condensed matter systems. 
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1 Introduction 


Conformal field theories (CFTs) lie at the heart of theoretical physics, describing critical 
phenomena in statistical and condensed matter systems, quantum gravity via the AdS/CFT 
correspondence, and possible solutions to the hierarchy problem (and other puzzles) in 
physics beyond the standard model. Quite generally, they serve as the endpoints of renor¬ 
malization group flows in quantum held theory. The conformal bootstrap [1, 2] aims to use 
general consistency conditions to map out and solve CFTs, even when they are strongly- 
coupled and do not have a useful Lagrangian description. 

In recent years great progress has been made in the conformal bootstrap in d > 2, 
including rigorous bounds on operator dimensions and operator product expansion (OPE) 
coefficients [3-32], analytical constraints [33-45], and methods for approximate direct solu¬ 
tions to the bootstrap [46-49], including a precise determination of the low-lying spectrum in 
the 3d Ising model under the conjecture that the conformal central charge is minimized [50]. 
These results have come almost exclusively from analyzing 4-point correlation functions of 
identical operators. It is tantalizing that even more powerful constraints may come from 
mixed correlators. 

In [51] some of the present authors demonstrated that semidehnite programming tech¬ 
niques can very generally be applied to systems of mixed correlators. In 3d CFTs with a Z 2 
symmetry, one relevant Z 2 -odd operator a, and one relevant Z 2 -even operator e, the mixed 
correlator bootstrap leads to a small and isolated allowed region in operator dimension 
space consistent with the known dimensions in the 3d Ising CFT. With the assistance of 
improved algorithms for high-precision semidehnite programming [52], this approach has 
culminated in the world’s most precise determinations of the leading operator dimensions 
(A^, AQ = (0.518151(6), 1.41264(6)) in the 3d Ising CFT. 

The immediate question is whether the same approach can be used to rigorously isolate 
and precisely determine spectra in the zoo of other known (and perhaps unknown) CFTs, 
particularly those with physical importance. In this work we focus on 3d CFTs with 0{N) 
global symmetry, previously studied using numerical bootstrap techniques in [15, 22]. We 
will show that the CFTs known as the 0{N) vector models can be similarly isolated using 
a system of mixed correlators containing the leading 0{N) vector and singlet s, assumed 
to be the only relevant operators in their symmetry representations. 

We focus on the physically most interesting cases N = 2,3,4 (e.g., see [53]) where the 
large-iV expansion fails. We do additional checks at iV = 20. A summary of the constraints 
on the leading scaling dimensions found in this work are shown in hgure 1. We also make 
precise determinations of the current central charge (JJ) oc Cj for N = 2,3. This coefficient 
is particularly interesting because it describes conductivity properties of materials in the 
vicinity of their critical point [54]. 

The 3d 0(2) model (or XY model) has a beautiful experimental realization in superfluid 
^He [55] which has yielded results for A^ that are in ~ Sa tension with the leading Monte 
Carlo and high temperature expansion computations [56]. Our results are not yet precise 
enough to resolve this discrepancy, but we are optimistic that the approach we outline in 
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this work will be able to do so in the near future. More generally, the results of this work 
give us hope that the same techniques can be used to to solve other interesting strongly- 
coupled CFTs, such as the 3d Gross-Neveu models, 3d Chern-Simons and gauge theories 
coupled to matter, 4d QCD in the conformal window, AA = 4 supersymmetric Yang-Mills 
theory, and more. 

The structure of this paper is as follows. In section 2, we summarize the crossing 
symmetry conditions arising from systems of correlators in 3d CFTs with 0{N) symmetry, 
and discuss how to study them with semidefinite programming. In section 3, we describe 
our results and in section 4 we discuss several directions for future work. Details of our 
implementation are given in appendix A. An exploration of the role of the leading symmetric 
tensor is given in appendix B. 
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Figure 1: Allowed regions for operator dimensions in 3d CFTs with an 0{N) global symmetry 
and exactly one relevant scalar (pi in the vector representation and one relevant scalar s in 
the singlet representation of 0{N), for N = 1,2,3,4,20. The case A = 1, corresponding to 
the 3d Ising model, is from [51]. The allowed regions for N = 2,3,4,20 were computed with 
A = 35, where A (defined in appendix A) is related to the number of derivatives of the crossing 
equation used. Each region is roughly triangular, with an upper-left vertex that corresponds 
to the kinks in previous bounds [15]. Further allowed regions may exist outside the range of 
this plot; we leave their exploration to future work. 
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2 Crossing Symmetry with Multiple Correlators 


Let us begin by summarizing the general form of the crossing relation for a collection of 
scalar helds (pi = (01, 02, 03, • • •)• We take the 0j to have dimensions Aj and for the moment 
we do not assume any symmetry relating them. Taking the OPE of the hrst two and last 
two operators, the 4-point function looks like: 


(0i(a:i)0j(a;2)0A:(a:3)0z(x4)) = 

1 


/ Xi4^ 

^kl 

) XijoXkiog^p 

0 

(m,u). 

Ai+Aj Ak+Ai 
Xi2 -('34 

Vaq4/ 

Ul3/ 

u = 

9 9 5 ^ 

,^2 2 
X 14 JP 23 

9 9 5 



( 2 , 1 ) 


X 13 X 24 

X 13 X 24 





where Xij = \xi — Xj\, Aij = Aj — A^, and u, v are the standard conformal invariants. 
The subscripts A,i refer to the dimension and spin of the operator O. We refer to [51] for 
details about how to compute the conformal blocks in any dimension and for 

arbitrary values of Aij. We also have the symmetry property Ajjo = {—lYXjio- 

Crossing symmetry of the correlation function requires that OPEs taken in different 
orders must produce the same result. As an example, exchanging (l,i) -H- (3, fc) gives the 
conditions: 


V 2 


o 


A^+Aj 

= U 2 


o 


( 2 . 2 ) 


It is convenient to symmetrize/anti-symmetrize in m,u, which leads to the two equations: 


0 = \^\ijoXkioFY]^Yu,v) ± XkjoXiioFY’f Y'^.v) 


(2.3) 


where 


o 


Up,A,((“.'") =« ” Sa,i (“,<') =F“ " 9a, i (v,u). 


(2.4) 


The functions are symmetric under exchanging i k and j 1. 


2.1 0{N) Models 

We now restrict our discussion to the case where 0* transforms in the vector representation 
of a global 0{N) symmetry. When the helds entering the four-point function are charged 
under global symmetries, the conformal block expansion can be organized in symmetry 
structures corresponding to irreducible representations appearing in the OPE 0j x 0^. This 
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gives the equations^ 


0 {5ij5ki i 5jk5ii) ^ ^ 

Os,i+ 


2 Tp4’<t>,4’4 


+ 



OT,e+ 



OA,e- 

which lead to three independent snm rnles after reading off the coefficients of each index 


strnctnre. Here, Os, Ot, Oa denote operators in the singlet, traceless symmetric tensor, 
and antisymmetric tensor representations of 0{N), O' refers to operators with even spin, 
and O refers to odd spin. The snm over spins is determined by the symmetry properties of 
the representations nnder exchange of two indices. 

In what follows, we will nse s, s', s" ,... to refer to the singlet scalars in increasing order 
of dimension. For example, s is the lowest-dimension singlet scalar in the theory. Similarly, 
and (j),(j)' ,(j)",... refer to scalars in the traceless symmetric tensor and vector 
representations, in increasing order of dimension. 

We would like to supplement the above equations with crossing symmetry constraints 
from other four-point functions. The simplest choice is to consider all nonvanishing four- 
point functions of 0* with the lowest dimension singlet scalar operator s. Another interesting 
choice would be the lowest dimension scalar in the traceless symmetric tensor representation 
tij. However the OPEs tij x tki and Uj x 0^ contain many additional 0{N) representations, 
increasing the complexity of the crossing equations. We leave the analysis of external tij 
operators to the future. 

Thus we consider the four-point functions {(pOjSs) and (ssss), which give rise to four 
additional sum rules after grouping the terms with the same index structure. In total this 

^Note that we are following the conformal block conventions of [51], which contain a factor of (—1)^ 
relative to the conventions used in the previous global symmetry studies [10, 15]. This leads to a different 
sign in front of the contributions of the Oa operators. 
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leads to a system of seven equations: 
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Note that the final line represents two equations, corresponding to the choice of ±. We can 
rewrite these equations in vector notation as 


0 


( ^<f>4>Os ^ssOs ) Vs,A,e 

Os/+ 


\<f>Os 

^ssOs 


+ 


E 

OT,e+ 




+ 

Oa/- Ov/^ 


(2.7) 


where Vr, Va, Vv are a 7-dimensional vectors and Vs is a 7-vector of 2 x 2 matrices: 
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( 2 . 8 ) 
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2.1.1 A Note on Symmetries 

We are primarily interested in theories with 0{N) symmetry. However, our bounds will 
also apply to theories with the weaker condition of SO{N) symmetry. This point deserves 
discussion. 

The group 0{N) includes reflections, so its representation theory is slightly different 
from that of SO{N). In particular is not an invariant tensor of 0{N) because it 

changes sign under reflections. For odd N = 2k + 1, 0{2k + 1) symmetry is equivalent to 
SO(2k + 1) symmetry plus an additional Z 2 symmetry. For even N = 2k, the orthogonal 
group is a semidirect product 0(2k) = Z 2 ix S0(2k), so it is not equivalent to an extra Z 2 . 

Let us consider whether the crossing equations must be modihed in the case of only 
SO{N) symmetry. In theories with SO{2) symmetry, the antisymmetric tensor represen¬ 
tation is isomorphic to the singlet representation. (This is not true for 0(2) because the 
isomorphism involves eij.) However in the crossing equation (2.7), antisymmetric tensors 
appear with odd spin, while singlets appear with even spin. Thus, the coincidence between 
A and S does not lead to additional relations in (2.7). 

For theories with SO (3) symmetry, the antisymmetric tensor representation is equivalent 
to the vector representation. Thus, antisymmetric odd spin operators appearing in 0 x 0 
may also appear in 0 x s. This does not affect (2.7) because there is no a priori relationship 
between and X(f,sO- However, it is now possible to have a nonvanishing four-point 

function {(f)i(f)j(f)ks) proportional to e^jk- Including crossing symmetry of this four-point 
function cannot change the resulting dimension bounds without additional assumptions. 
The reason is as follows. Any bound computed from (2.7) without using crossing of (000s) 
is still valid. Hence, the bounds cannot weaken. However, because any 0(3)-invariant theory 
is also S'0(3)-invariant, any bound computed while demanding crossing of (000s) must also 
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apply to 0(3)-invariant theories. So the bounds cannot strengthen. Crossing for {(fxjxps) 
only becomes important if we input that is nonzero for a particular operator.^ 

This would guarantee our theory does not have 0(3) symmetry. 

For 5*0 (4), the new ingredient is that the antisymmetric tensor representation can be 
decomposed into self-dual and anti-self-dual two-forms. As explained in [10], this leads to 
an additional independent sum rule 

E E = 0. 

A+,^- A-,i- 

where A± represent self-dual and anti-self-dual operators. By the same reasoning as in 
the case of S'0(3), this crossing equation cannot affect the bounds from (2.7) without 
additional assumptions. We can also see this directly from (2.10) together with (2.7): in 
the semidehnite program used to derive operator dimension bounds, we may always take 
the functional acting on (2.10) to be zero. An exception occurs if we know an operator is 
present with XtjupoA 0 but — 0 versa). Then we can include that operator 

with other operators whose OPE coefficients are known (usually just the unit operator) and 
the resulting semidehnite program will be different. 

For SO{N) with A^ > 5, no coincidences occur in the representation ring that would 
be relevant for the system of correlators considered here. In conclusion, (2.7) and the 
semidehnite program discussed below remain valid in the case of SO{N) symmetry. Bounds 
on theories with SO{N) symmetry can diher only if we input additional information into 
the crossing equations that distinguishes them from 0(A^)-invariant theories (for example, 
known nonzero OPE coefficients). 

2.2 Bounds from Semidefinite Programming 

As explained in [51], solutions to vector equations of the form (2.7) can be constrained using 
semidehnite programming (SDP). We refer to [51] for details. Here we simply present the 
problem we must solve. To rule out a hypothetical CFT spectrum, we must hnd a vector 
of linear functionals a = (ai, 02 ,..., aj) such that 


1 )«• ^5,0,0 ( J ^ 

1 > 0, for the identity operator. 

(2,11) 

•3 ■ f^r,A,£ ^ 0, 

for all traceless symetric tensors with i even. 

(2.12) 

o' 

Al 

<f 

for all antisymmetric tensors with i odd. 

(2.13) 

0 • Vv,A,e > 0, 

for all 0{N) vectors with any i, 

(2.14) 

« • Vs, A/ — 0) 

for all singlets with £ even. 

(2.15) 


Here, the notation “E 0” means “is positive semidehnite.” If such a functional exists for a 
hypothetical CFT spectrum, then that spectrum is inconsistent with crossing symmetry. In 

^In practice, this means we would group this operator with the unit operator and other operators whose 
OPE coefficients are known in the semidefinite program. 
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addition to any explicit assnmptions placed on the allowed valnes of A, we impose that all 
operators must satisfy the unitarity bound 



( 2 . 16 ) 


A 


where D = 3 is the spacetime dimension. 

Additional information about the spectrum can weaken the above constraints, making 
the search for the functional a easier, and further restricting the allowed theories. A few 
specihc assumptions will be important in what follows: 

• The 3d 0{N) vector models, which are our main focus, are believed to have exactly 

one relevant singlet scalar s, 0{N) vector scalar 0,, and traceless symmetric scalar 
tij.^ We will often assume gaps to the second-lowest dimension operators in 

each of these sectors. These assumptions affect (2.12), (2.14), and (2.15). 

• Another important input is the equality of the OPE coefficients = X^g^. This 
is a trivial consequence of conformal invariance. It is important that 0 and s be 
isolated in the operator spectrum for us to be able to exploit this constraint. For 
instance, imagine there were two singlet scalars si ,2 with the same dimension. Then 

"would appear in (2.7). This combination does not satisfy 

\fake \ 

• We will sometimes assume additional gaps to derive lower bounds on OPE coefficients. 
For instance, to obtain a lower bound on the coefficient of the conserved 0{N) current 
in the 0, x (pj OPE, we will need to assume a gap between the hrst and second spin-1 
antisymmetric tensor operators. 

As an example, (2.17) shows a semidehnite program that incorporates symmetry of X^^g 
and the assumption that 0*, s are the only relevant scalars in their respective sectors: 



(unit operator) 


<3 • VA,A,e A 0, 

<3 • Vv,A/ — 


CK • Vt,a,(. a 0, 


A > £ = 0, and 

/S.> ^ + D — 2, £>0 even; 
^> ^ + D — 2, i odd; 


A > 
A > 
A > 


i + D-2, 
i + D-2, 


D-2 
2 ’ 


A> D £ = 0, and (2-17) 

A>£ + D-2, £>0; 


3 • Vs,A,i ^ 0 ) 


A > D, £ = 0, and 

A > £ + D — 2 £ > 0 even; 



^Additional relevant scalars could be present in other representations. 
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The final constraint in (2.17) imposes the appearance of (j)i, s in the OPEs and incorporates 
the eqnality = X^s<i>-'^ It replaces two otherwise independent constraints on Vs and 
Vy. As previously mentioned, if we assume no gap between (pi, s and the next operators in 
each sector, enforcing symmetry of the OPE coefficients will have no effect: indeed each of 
the terms in this constraint would be independently positive-semidefinite, since the other 
inequalities imply a • q ^ 0 and a ■ 17y,A^+5,o > 0 for 5 arbitrary small. 

Finally, one might want to enforce the existence of a unique relevant scalar operator, 
with dimension A^, transforming in the traceless symmetric representation. In this case the 
symmetric tensor constraint is replaced by 

a • Vt,a/ — 0) A = At or A > D, £ = 0, and 

A>i + D-2, i>0even. (2.18) 


3 Results 

3.1 0(2) 

To begin, let us recall the bounds on A^, Ag computed in [15] using the correlation function 
{pipjpk4>i) (see figure 2). Like the Ising model bounds computed in [12, 50], this single¬ 
correlator bound has an excluded upper region, an allowed lower region, and a kink in the 
curve separating the two. The position of this kink corresponds closely to where we expect 
the 0(2) model to lie, and one of our goals is to prove using the bootstrap that the 0(2) 
model does indeed live at the kink.® If we assume that s is the only relevant 0(2) singlet, 
then a small portion of the allowed region below the kink gets carved away, analogous to 
the Ising case in [51]. 

Adding the constraints of crossing symmetry and unitarity for the full system of correla¬ 
tors {pppp), {ppss), {ssss) does not change these bounds without additional assumptions. 
However, having access to the correlator (ppss) lets us input information special to the 0(2) 
model that does have an effect. We expect that p is the only relevant 0(2) vector in the 
theory. One way to understand this fact is via the equation of motion at the Wilson-Fisher 
fixed point in 4 — e dimensions. 


npi oc Xp'^pi- (3.1) 

This equation implies that the operator is a descendent, so there is a gap in the spectrum 
of 0(2)-vector primaries between pi and the next operator in this sector, which is a linear 
combination of pip'^ and ppdpY. The equation of motion makes sense in perturbation 
theory e <C 1. However, it is reasonable to expect gaps in the spectrum to be robust as 

^In writing this constraint, we have assumed the scalar conformal blocks are normalized so that 
to leading order in u, where C is a A-independent constant. 

®The sharpness of the kink depends on the number of derivatives A used when computing the bound 
(appendix A). Figure 2 was computed at a lower derivative order than we use for most of this work, so the 
kink is relatively smooth. 
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Figure 2: Allowed region for (A^,As) in 3d CFTs with 0(2) symmetry. The light blue 
region makes no additional assumptions and was computed in [15] using the correlator {(f)(l>(f)(l>) 
at A = 19. The medium blue region was computed from the system of correlators {(f>(p(j)(p), 
{(jxpss), (ssss) at A = 19, and assumes A^ and A^ are the only relevant dimensions in the 
vector and singlet scalar channels at which contributions appear. The dark blue region is 
computed similarly, but additionally assumes the OPE coefficient relation X(j)(j)s = This 

latter assumption leads to a small closed region in the vicinity of the red cross, which represents 
the Monte Carlo estimate for the position of the 0(2) model from [56]. 
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e gets larger. In particular, we expect this gap to persist as e —)■ 1. Thus, a gap in the 
0(2)-vector sector reflects the equations of motion of the 0(2) model. 

We do not know if there is sharp experimental evidence for the claim that the 0(2) model 
contains exactly one relevant 0(2)-vector scalar. The cleanest experimental realization of 
the 0(2) model is the superfluid transition in '^He [55]. This theory has microscopic 0(2) 
symmetry, so one cannot easily determine the number of relevant 0(2)-vector scalars by 
counting order parameters. The number could be determined by counting order parameters 
in systems where the 0(2) symmetry is emergent. 

As explained above, it is natural to impose a gap in both the 0(2) vector and singlet 
sectors in our formalism, giving rise to the medium blue region in figure 2. Another 
important constraint is symmetry of the OPE coefficient Adding this constraint 

gives the dark blue region in figure 2; a close-up view of the 0(2) model point is shown in 
figure 3, which we show for increasing numbers of derivatives A = 19, 27, 35 (see appendix A). 
We now have a closed island around the expected position of the 0(2) model, very close to 
the original kink in figure 2. The bounds strengthen as A increases. However, the allowed 
regions apparently do not shrink as quickly as in the case of the 3d Ising GET [52]. Thus, our 
determination of (A^, A^) is unfortunately not competitive with the best available Monte 
Carlo [56] and experimental [55] results (though it is consistent with both).® We conjecture 
that including additional crossing relations (such as those involving the symmetric tensor 
tij) will give even stronger bounds; we plan to explore this possibility in future work. 

In addition to gaps in the 0(2)-vector and singlet sectors, we also expect that the 
0(2) model has a single relevant traceless symmetric tensor tij. Let us finally impose this 
condition by demanding that fC has dimension above D = 3 and scanning over A^ along with 
A^, As. The result is a three-dimensional island for the relevant scalar operator dimensions, 
which we show in figure 4. Our errors for the symmetric-tensor dimension A^ are much more 
competitive with previous determinations. By scanning over different values of (A^, A^) in 
the allowed region and computing the allowed range of A^ at A = 35, we estimate 

1.2325 < At < 1.239 (0(2) model), (3.2) 

which is consistent with previous results from the pseudo-e expansion approach [57] giving 
At = 1.237(4). 

3.2 0{N),N>2 

The bounds for N > 2 are similar to the case of A = 2. In figure 5, we show the allowed 
region of (A,^, A^) for theories with 0(3) symmetry, assuming 0 and s are the only relevant 
scalars in their respective 0(N) representations, and using symmetry of the OPE coefficient 
A^^s- We expect that an additional scan over At would yield a 3d island similar to figure 4. 

®Note that '’’He experiments cannot easily determine because the 0(2) symmetry is realized 
microscopically. Some results constraining have been reported from NMR experiments (e.g., as 
summarized in [53]) but they are not very precise. 
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0(2): Multi-Correlator with OPE relation 
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Figure 3: Allowed regions for (A^jA^) in 3d CFTs with 0(2} symmetry and exactly one 
relevant 0(2) vector (j) and singlet s, computed from the system of correlators {((xpc^xj)), {(fxpss), 
and (ssss) using SDPB with A = 19, 27, and 35 (see appendix A). The smallest region (darkest 
blue) corresponds to A = 35. The green rectangle represents the Monte Carlo estimate 
[56]. The red lines represent the Icr (solid) and 3cr (dashed) confidence intervals for A^ from 
experiment [55]. The allowed/disallowed regions in this work were computed by scanning over 
a lattice of points in operator dimension space. For visual simplicity, we fit the boundaries 
with curves and show the resulting curves. Consequently, the actual position of the boundary 
between allowed and disallowed is subject to some error (small compared to size of the regions 
themselves). We tabulate this error in appendix A. 
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Multi-Correlator 3D Scan 



Figure 4: Allowed region (orange) for (A^,As,At) in a 3d CFT with 0(2) symmetry and 
exactly one relevant 0(2)-vector cf), 0(2) singlet s, and 0(2) traceless symmetric-tensor t. This 
region was computed using SDPB with A = 19. The green rectangle represents the error bars 
from Monte Carlo [56] and the pseudo-e expansion approach [57]. Note that our estimate for 
At in (3.2) was computed with A = 35, so it is more precise than the region pictured here. 
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Figure 5: Allowed regions for (A^,As) in 3d CFTs with 0(3) symmetry and exactly one 
relevant 0(3)-vector cj) and 0(3) singlet s, computed using SDPB with A = 19,27, and 35 (see 
appendix A). The smallest region (darkest blue) corresponds to A = 35. The green rectangle 
represents the Monte Carlo estimate [58]. 


By performing this scan at a few values of (A^, A^), we estimate 


((9(3) model), 


(3.3) 


1.204 < A, < 1.215 


which is consistent with previous results from the pseudo-e expansion approach [57] giving 


At = 1.211(3). 


In hgure 6, we show the allowed region of (A^, A^) for the (9(4) model, with the same 
assumptions as discussed above for (9(3). A clear trend is that the allowed region is growing 
with N. For example, at A = 19, the (9(4) allowed region isn’t even an island — it connects 
to a larger region not shown in the plot. Increasing the number of derivatives to A = 35 
shrinks the region, but it is not as small as in the case of (9(2) or (9(3). 

The trend of lower-precision determinations at larger N reverses at some point. For 
example, in figure 1, the allowed region for A = 20 is smaller again than the (9(4) region. 
The relative size of the (9(4) region and the (9(20) region is A-dependent, and we have not 
studied the pattern for general N in detail. However, as an important check we note that 
the (9(20) island in figure 1 is nicely compatible with the 1/N expansion (see [15]), giving 
the point (A<^, A^) ~ (.5064,1.938) which sits in the upper-left corner of the allowed region. 

Finally, we remark that all of the constraints on operator dimensions found above can 
be reinterpreted in terms of constraints on critical exponents. Following standard critical 
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Figure 6: Allowed regions for (A^,As) in 3d CFTs with 0(4) symmetry and exactly one 
relevant 0(4)-vector cj) and 0(4) singlet s, computed using SDPB with A = 19,27, and 35 (see 
appendix A). The smallest region (darkest blue) corresponds to A = 35. The green rectangle 
represents the Monte Carlo estimate [59]. 


exponent notation (see [53]), the relations are given by 


T] = 2A^ - 1 , 


/S 


Aip 

3-A,’ 


1 ^ _ 3 - 2A0 

“ 3-A,’ ^ “ 3-A, 

_ 3 - A^ 3 - 4A^ 

A^ ’ ^ 3-A, 


a = 
02 = 


3-2A, 
3-A, ’ 
3-Ai 
3-A/ 


(3.4) 


3.3 Current Central Charges 

Let be the conserved currents that generate 0{N) transformations. Jj^ (x) is an 0{N) 

antisymmetric tensor with spin 1 and dimension 2. Its 2-point function is determined by 
conformal and 0{N) symmetry to be 






c, 


(dvr 


\2 ^4 


X 


12 L 


V 


uu X2Y{Xi- X2y 


•^12 


(3.5) 


We call the normalization coefficient Cj from Eq. (3.5) the current central charge.^ The con¬ 
served current J^j appears in the sum over antisymmetric-tensor operators Oa in Eq. (2.7). 

^This name is by analogy with the case of 2d CFTs, where Cj appears as a central element in an affine 
Kac-Moody algebra. In higher dimensional CFTs, Cj is not an element of a nontrivial algebra in general, 
though it can be in special cases [37]. 
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A Ward identity relates the OPE coefficient Xjcj,^ to Cj. In our conventions 


’ (3-6) 

where = 2 is the free theory value of Cj [60, 61]. In the 0{N) vector models Cj is 
known to have the large N and e expansions [62] 


Cj 


^free 


(i=3 


32 1 

^iV 


+ 0 



Cj 


^free 


d=A—e 


1 - 


3(iV + 2) 2 

4(iV + 8)2^ 


+ 0 ( 6 =^). 


(3.7) 


Note that both of these expansions predict that Cj will be smaller than the free value. 


It is well known that the conformal bootstrap allows one to place upper bounds on OPE 
coefficients, or equivalently a lower bound on Cj. Previously such bounds were explored in 
d = 4 in [6, 10] and in d = 3, 5 in [22]. To hnd such a bound, we search for a functional a 
with the following properties (cf. eq. (2.17)): 


d ■ 


— 1, 

(normalization) 

Oi ■ Vt,A,^ — 

A > ^-2 


A>7 + T)-2, 

d ■ Vaa.i a 0, 

A>7 + T)-2, 

<7 • ^y,A,£ > 0, 

A>D 


A>7 + T)-2, 

<7 • ^s,A,£ ^ 0, 

A>D, 


A> l + D-2 




i = 0, and 
i > 0 even; 
i odd; 

7 = 0, and 

7 > 0 ; 

7 = 0, and 
7 > 0 even; 


(3.8) 


Notice that compared to (2.17), we have dropped the assumption of the functional a being 
positive on the identity operator contribution and we chose a convenient normalization for 
a. It follows then from the crossing equation (2.7) that 

^^<-(1 l)S-ho,o(]). (3.9) 

Therefore, Ending a functional a sets a lower bound on Cj. To improve the bound, we 
should minimize the RHS of (3.9). We thus seek to minimize 

-(11 )«-E5,o,o( J ) , (3.10) 

subject to the constraints (3.8). This type of problem can be efficiently solved using SDPB. 
In this way, we set a lower bound on Cj for all allowed values of A^, Ag. 

We can also set an upper bound on Cj, provided we additionally assume a gap in the 
spin-1 antisymmetric tensor sector. At this point it is not clear what gap we should assume, 
but to stay in the spirit of our previous assumptions, we will assume that the dimension of 
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the second spin-1 antisymmetric tensor satisfies Aj/ > 3, so that the current J^ - is the only 
relevant operator in this sector. We now search for a functional a (different from the one 
above) that satisfies 


a-14,A,i>0, A >3, (3.11) 

a-14,A,^>0, A >£ + /!-2, £>lodd, (3.12) 


and is normalized so that 


= (3.13) 

The constraints on a coming from the singlet and traceless symmetric-tensor sectors stay 
the same as in (3.8). An upper bound on Cj then follows from (2.7): 

l)3'W.(;). (3.14) 

Our upper and lower bounds on Cj, expressed as a function of A^ and A^, are shown 
in figures 7 and 8 for 0(2) and 0(3) symmetry, respectively. The allowed region for a given 
N consists of a 3d island in (A^, A^, Oj) space. This determines the current central charge 
to within the height of the island. For the two physically most interesting cases, N = 2 and 
A = 3, we find: 

N = 2: ^ = 0.9050(16), A = 3: ^ = 0.9065(27). (3.15) 

As an additional check, we also computed Cj for N = 20: 

A = 20 : ^ = 0.9674(8). (3.16) 

This result agrees within 0.5% accuracy with the leading 1/A expansion result, Cj/Cj*^^ ~ 

0.964 [62]. 

Recently, the current central charge attracted some interest in studies of transport 
properties of 0(A) symmetric systems near a quantum critical point, where Cj can be 
related to the conductivity at zero temperature. In particular, using the OPE it was found 
in [54] that the asymptotic behavior of conductivity at low temperature is given by 

+ ^ -i2ACtiH,,(^ +..., (3.17) 

where aq = /h is the conductance quantum. Here, (Too is the (unitless) conductivity at 
high frequency and zero temperature which is related to Oj as 

cToo = Oj/32. (3.18) 
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0(2): Cj allowed region 


0(2): Projected Cj allowed region 



A 


Figure 7: The left panel shows the allowed values of Cj as a function of and in 0(2) 
symmetric theories. The right panel is the projection of the allowed region onto the (A^, Cj) 
plane. Both plots are computed using SDPB with A = 27. 


0(3): Cj allowed region 


0(3): Projected Cj allowed region 
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Figure 8: The left panel shows the allowed values of Cj as a function of A^ and A^ in 0(3) 
symmetric theories. The right panel is the projection of the allowed region onto the (A^, Cj) 
plane. Both plots are computed using SDPB with A = 27. 









Furthermore, Ct is the central charge of the theory, C is the (JJs) OPE coefficient, and 7 
is one of the (JJT) OPE coefficients, where T is the energy-momentum tensor. B and H^x 
are the hnite temperature one-point function coefficients: 

(s)t = , {Txx:)t = Hxx^T^. (3.19) 

Of all the GET data that goes into (3.17), we have determined (Too and for the 0{N) 
vector models in this work, while Ct was estimated using bootstrap methods before in [15]. 
The OPE coefficients C and 7 can not be determined in our setup, but could in principle 
be obtained by including the conserved current as an external operator in the crossing 
equations. The one-point functions B and Hxx are in principle determined by the spectrum 
and OPE coefficients of the theory [63]. However, to compute them we would need to 
know the high-dimension operator spectrum. This is still out of the reach of the conformal 
bootstrap approach. 

Of particular interest for physical applications is the N = 2 case, which describes 
superfluid-insulator transitions in systems with two spatial dimensions [64, 65]. Some 
examples of such systems are thin Elms of superconducting materials, Josephson junction 
arrays, and cold atoms trapped in an optical lattice. In these systems the parameter aoo 
is the high-frequency limit of the conductivity. This quantity has not yet been measured 
in experiments, but was recently compnted in Quantum Monte Carlo simulations [66, 54], 
[67], and [68] to be 27r(T^'^ = 0.3605(3), 0.359(4), and 0.355(5), respectively®. Our rigorous 
result 27r(T^°°*®*’’®'P = 0.3554(6) is in excellent agreement with these determinations and is 
significantly more precise after systematic uncertainties are taken into acconnt. 


4 Conclusions 

In this work, we used the conformal bootstrap with multiple correlators to set more stringent 
bounds on the operator spectrum of 3d CFTs with 0{N) symmetry. The multiple correlator 
approach works in this setting similarly to the case of Z 2 -symmetric CFTs - including 
mixed correlators opens access to parts of the spectrum that are inaccessible with a single 
correlator. In this work we considered mixed correlators of an 0{N) singlet and an 0{N) 
vector, gaining access to the sector of 0{N) vectors. We can then additionally input 
assumptions about the operator spectrum in that sector. As a result, we exclude large 
portions of the allowed space of CFTs. This reaffirms conclusions from previous works on 
the 3d Ising model; it is important and fruitful to consider multiple crossing equations. We 
believe that including mixed correlators will be rewarding in many other bootstrap studies 
that are currently ongoing. 

Specifically, for 0{N) symmetric CFTs, we found that the scaling dimensions of the 
lowest 0{N) vector scalar 0 and 0{N) singlet scalar s are constrained to lie in a closed 
region in the (A(^,As) plane. Our assumptions, besides conformal and 0{N) symmetry, 

®These uncertainties reflect statistical errors but may not include systematic effects, conservatively 
estimated in [66] to be 5 — 10%. We thank Subir Sachdev, Erik S0rensen, and William Witczak-Krempa for 
correspondence on this point. 
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were crossing symmetry, unitarity, and — crucially — the absence of other relevant scalars 
in the 0{N) singlet and vector sectors. This is completely analogous to the Z 2 -symmetric 
case where similar assumptions isolate a small allowed region around the Ising model in 
the (Ao-jAg) plane. Our allowed regions represent rigorous upper and lower bounds on 
dimensions in the 0{N) models. In principle, this approach could be used to compute the 
scaling dimensions of the 0{N) models to a very high precision, assuming that the allowed 
region will shrink to a point with increased computational power. However, our results 
suggest that the region either does not shrink to a point, or the convergence is slow in the 
present setup. Therefore, our uncertainties are currently larger than the error bars obtained 
using other methods.® In particular, we have not yet resolved the discrepancy between 
Monte Carlo simulations and experiment for the value of A,, in the 0(2) model. 

Including more correlators could result in signihcantly improved bounds on operator 
dimensions. In the case of 0{N) symmetric CFTs, it would be natural to include the lowest 
dimension 0{N) symmetric tensor as an external operator in the crossing equations. In the 
0{N) models, this operator actually has a lower dimension than s. This is an important 
difference from the Ising model, where 0 and s are the two lowest dimensional scalars in any 
sector of the theory. Our present bounds on the lowest symmetric tensor treated it as an 
internal operator in the crossing equations. Including it as an external operator would open 
access to many other 0{N) representations. Perhaps the 0{N) models are not uniquely 
determined by the condition of only one relevant 0{N) singlet and vector scalar, and we 
must also specify something about these other representations. Studying the 0{N) models 
in other dimensions (such as in 5d [22, 28, 30, 69, 70]) may also help to shed light on these 
issues. We plan to further explore these questions in the future. 

In addition to scaling dimensions, it is also important to determine OPE coefficients. 
Here we presented an example in the computation of the current central charge Cj. In 
the case of 0(2) symmetry, this yields the current most precise prediction for the high- 
frequency conductivity in 0(2)-symmetric systems at criticality. It will be interesting to 
extend these mixed-correlator computations to other OPE coefficients in the 0{N) models 
such as the stress-tensor central charge Ct and 3-point coefficients appearing in (JJs) and 
(JJT) which control frequency-dependent corrections to conductivity. Pursuing the latter 
will require implementing the bootstrap for current 4-point functions, a technical challenge 
for which efforts are ongoing in the bootstrap community. 

More generally, the results of this work make it seem likely that scaling dimensions 
in many other strongly-interacting CFTs can be rigorously determined using the multiple 
correlator bootstrap. It will be interesting to study mixed correlators in 3d CFTs with 
fermions and gauge helds - it is plausible that similar islands can be found for the 3d Gross- 
Neveu models and 3d Chern-Simons and gauge theories coupled to matter. In 4d, we hope 
that by pursuing the mixed correlator bootstrap we will eventually be able to isolate and 
rigorously determine the conformal window of QCD. It also be interesting to apply this 
approach to theories with conformal manifolds to see the emergence of lines and surfaces 
of allowed dimensions; a concrete application would be to extend the analysis of [14, 23] to 

®If one is willing to assume that the 0{N) models live near kinks in our allowed regions, then more 
precise determinations are possible. 


22 



mixed correlators and pursue a rigorous study of the dimension of the Konishi operator in 
AA = 4 supersymmetric Yang-Mills theory at hnite N. The time is ripe to set sail away from 
our archipelago and explore the vast ocean of CFTs! 
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A Implementation Details 

As described in [51], the problem of Ending a satisfying (2.11) can be transformed into 
a semidehnite program. Firstly, we must approximate derivatives of Vs, Ur, U 4 , and Vy as 
positive functions times polynomials in A. We do this by computing rational approximations 
for conformal blocks using the recursion relation described in [51]. Keeping only the 
polynomial numerator in these rational approximations, (2.11) becomes a “polynomial 
matrix program” (PMP), which can be solved with SDPB [52]. 

Three choices must be made to compute the PMP. Firstly, k (dehned in appendix A 
of [52]) determines how many poles to include in the rational approximation for conformal 
blocks. Secondly, A determines which derivatives of conformal blocks to include in the 
functionals a. Specifically, we take 


a.(f) = 5^ (A.l) 

m+n<A 

Some of these derivatives vanish by symmetry properties of F. The total number of nonzero 
components of a is 

+ , .L^J(L^J + i) 

^ 2 + 2 

Finally, we must choose which spins to include in the PMP. The number of spins depends 
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on A as follows 

^a=i9 = {0,..., 26} U {49, 50}, 

= {0,..., 26} U {29, 30, 33, 34, 37, 38,41,42,45,46,49, 50}, 

^^^35 = {0,..., 44} U {47,48, 51, 52, 55, 56, 59, 60, 63, 64, 67, 68}, 

= {0,..., 54} U {57, 58, 61, 62, 65, 66, 69, 70, 73, 74, 77, 78}. (A.3) 

We use Mathematica to compute and store tables of derivatives of conformal blocks. 
Another Mathematica program reads these tables, computes the polynomial matrices cor¬ 
responding to the Id’s, and uses the package SDPB. m to write the associated PMP to an xml 
hie. This xml hie is then used as input to SDPB. Our settings for SDPB are given in table 1. 


A 

19 

27 

35 

39 

K 

14 

20 

30 

36 

spins 

<S'a=19 

*S'a=27 

5'a=35 

*S'a=39 

precision 

448 

576 

768 

896 

findPrimalFeasible 

True 

True 

True 

True 

findDualFeasible 

True 

True 

True 

True 

detectPrimalFeasibleJump 

True 

True 

True 

True 

detectDualFeasibleJump 

True 

True 

True 

True 

dualityGapThreshold 

10-30 

O 

1 

CO 

o 

o 

1 

CO 

o 

10-70 

primalErrorThreshold 

o 

1 

CO 

o 

10-30 

O 

1 

4i^ 

O 

10-70 

dualErrorThreshold 

O 

1 

CO 

o 

10-30 

o 

1 

4i^ 

O 

10-70 

initialMatrixScalePrimal 

1040 

1050 

1050 

1060 

initialMatrixScaleDual 

1040 

1050 

1050 

o 

Cb 

o 

feasibleCenteringParameter 

0.1 

0.1 

0.1 

0.1 

infeasibleCenteringParameter 

0.3 

0.3 

0.3 

0.3 

stepLengthReduction 

0.7 

0.7 

0.7 

0.7 

choleskyStabilizeThreshold 

O 

1 

O 

O 

1 

O 

lO-ioo 

10-120 

maxC omp1ement arity 

lOioo 

10130 

10160 

o 

00 

o 


Table 1: SDPB parameters for the computations of scaling dimension bounds in this work. 
For Cj bounds we need to set all of the Boolean parameters in the table to False. In addition 
to that, we used dualityGapThreshold = 10“^*^, while all the rest of the parameters were 
kept at the same values as for the dimension bounds. 

Finally let us conclude with some comments on the precision of the plots presented in 
the main text. Conformal blocks of correlation functions involving operators of nonequal 
dimensions depend nontrivially on the diherence of the dimensions. Hence, when computing 
the boundary of various allowed regions, it is convenient to perform a scan over a lattice of 
points. The vectors generating the lattice points are shown in table 2. The smooth regions 
shown in hgs. 1, 3, 5, and 6 are the results of a least-squares £t, subject to the constraint 
that allowed lattice points should lie inside the curves while excluded ones lie outside. In 
table 2 we also show the maximal perpendicular distance of these points to the curves. 

The bounds on Cj shown in hgures 7 and 8 were computed for the lattices of points that 
were found to be allowed in hgures 3 and 5. For each point on the lattice, the bound on Cj 
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was determined to a precision of 10“^*^. The smooth regions were obtained by interpolation 
and the maximum distance of the computed points to the boundry of the shaded region is 
again reported in table 2. 




allowed 

excluded 

■^1 

V2 


A = 19 

0.00025 

0.00060 

(nr^yiir^) 

(0,10-^) 

0(2) 

A = 27 

0.000084 

0.00025 

(lo-^lo-^) 

(0,4-10-^) 


A = 35 

0.00021 

0.00062 

(5- 10-^5■ 10“^) 

o 

1 


A = 19 

0.00043 

0.0020 

O 

1 

o 

1 

(0,2-10-^) 

0(3) 

A = 27 

0.00044 

0.0019 

o 

1 

o 

1 

_ 

(0,2-10-3) 


A = 35 

0.00041 

0.0013 

o 

1 

o 

1 

(0,10-3) 


A = 19 

0.00040 

0.00041 

o 

1 

o 

1 

(0,2-10-3) 

0(4) 

A = 27 

0.00048 

0.00048 

o 

1 

o 

1 

(0,2-10-3) 


A = 35 

0.00029 

0.00062 

(10-^10-^) 

(0,2-10-3) 

0(20) 

A = 35 

0.00014 

0.00023 

o 

1 

o 

1 

4 ^ 

(0,2-10-3) 

0(2): Cj 

A = 27 

0.00005 

- 

o 

1 

o 

1 

(0,-10-3) 

0(3): Cj 

A = 27 

0.0001 

- 

o 

1 

o 

1 

to 

o 

1 

o: 


Table 2: Maximal distance between the computed allowed and excluded points and the curves 
shown in figs. 1, 3, 5, 6, 7 and 8. The vectors vi and V 2 describe the direction and spacing 
of the computed grids in the (A^, A^) plane. For the Cj bounds we use the same lattices in 
the (A^, A^) plane. The reported maximal distance in the table is the vertical distance of the 
computed points to the regions shown in the right panels of figures 7 and 8. 
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B Symmetric Tensor Scan 


In this appendix we collect some detailed scans of the allowed region of (A^, A^, A*) space 
for 0{N) models with N = 2,3,4. The results for the 0(2) model are also presented as a 
3d plot in hgure 4. Here we show plots in the (A,^, A^) plane at hxed values of A*. The 
scans for 0(2), 0(3) and 0(4) are shown in hgs. 9, 10, and 11, respectively. Blue points 
represent the allowed region at A = 19. The light blue shaded area is the allowed region 
at A = 35, but without any assumptions in symmetric tensor sector; those are the same 
allowed regions shown in hgs. 3, 5, and 6. The hnal allowed regions with the assumptions 
on At are thus given by the intersections of the dark blue and light blue regions. 

Qualitatively the picture is the same for each value of N and we expect that the 
projections of the 3d plot into the (A<^, A^) plane will look similar for even higher values of 
N. In particular, the lowest allowed values of At are obtained at the lower left corner of the 
allowed region in the (A<^, A^) plane, while the greatest values are obtained at upper right 
corner of the allowed region. This allows us to hnd general bounds on At without doing a 
whole scan over the (A^, A^) plane; it is enough to hnd bounds on At at the corner points. 
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Figure 9: Allowed points in the (A,^, A^) plane for different values of At in 0(2) symmetric 
CFTs at A = 19 (dark blue). The light blue shows the allowed region at A = 35 without 
any assumptions on the symmetric tensor spectrum. The green rectangle is the Monte Carlo 
estimate [56]. 
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Figure 10: Allowed points in the (A,^, A*) plane for different values of A^ in 0(3) symmetric 
CFTs at A = 19 (dark blue). The light blue shows the allowed region at A = 35 without 
any assumptions on the symmetric tensor spectrum. The green rectangle is the Monte Carlo 
estimate [58]. 
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Figure 11: Allowed points in the (A^, A*) plane for different values of A* in 0(4) symmetric 
CFTs at A = 19 (dark blue). The light blue shows the allowed region at A = 35 without 
any assumptions on the symmetric tensor spectrum. The green rectangle is the Monte Carlo 
estimate [59]. 
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